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We study analytically the features of the Cauchy horizon (CH) singularity inside a spherically- 
symmetric charged black hole, nonlinearly perturbed by a self-gravitating massless scalar field. We 
derive exact expressions for the divergence rate of the blue-shift factors, namely the derivatives in 
the outgoing direction of the scalar field $ and the area coordinate r. Both derivatives are found 
to grow along the contracting CH exactly like 1/r. Our results are valid everywhere along the CH 
singularity, up to the point of full focusing. These exact analytic expressions are verified numerically. 
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I. INTRODUCTION 

In the last few years, the investigation of spinning and 
charged black holes led to a new picture of the space- 
time singularity inside such black holes. According to 
this new picture, the Cauchy horizon (CH) evolves into 
a curvature singularity, which has the following two re- 
markable features: (i) It is null (rather than spacelike), 
(ii) It is weak (in Tipler's terminology namely, the 
tidal deformation experienced by extended physical ob- 
jects is finite at the null singularity. In the case of a 
spinning black hole, the evidence for the occurrence of 
the null weak singularity has emerged from a systematic 
linear and nonlinear perturbation analysis For the 
toy model of a spherical charged black hole, the main 
features of the singularity at the inner horizon were first 
deduced analytically from simplified models based on null 
fluids and later confirmed numerically for a model 

with a self-gravitating scalar field . (See also the ap- 
proximate leading-order analysis in |§|]). In addition, the 
local existence and genericity of the null weak singularity 
were shown mathematically in Ref . , and more recently 
(in the framework of plane-symmetric spacetimes) in Ref. 
JHJ. The compatibility of a null weak singularity with 
the constraint equations was shown in Ref. fllj| . 

We shall consider here the spherically-symmetric 
model of a charged black hole nonlinearly perturbed by a 
sclf-gravitating, minimally-coupled, massless scalar field. 
Despite its relative simplicity (compared to the analogous 
model of a spinning black hole), no systematic analytic 
study of this model has been carried out so far. The 
goal of this paper is to present a simple analytic calcula- 
tion, which may be the first step towards such a thorough 
analytic study: We quantitatively analyze the evolution 
of the divergent blue-shift factors along the contracting 
CH. It is well known (from both theoretical considera- 
tions and numerical simulations) that the singularity at 
the CH is characterized by finite values of the scalar field 
$ and the area coordinate r. (The latter is also known to 
decrease monotonically with increasing afnne parameter 
along the CH, due to the outflux of energy-momentum 



carried by the scalar field.) However, the gradients of 
$ and r diverge at the CH. More specifically, let V be 
a "Kruskal-like" ingoing null coordinate (i.e. an ingoing 
null coordinate for which the double-null metric function 
g uv is finite and nonvanishing at the Cauchy horizon - see 
below). Then, r v and $ v diverge at the CH. In this pa- 
per we shall calculate the evolution of r v and <E> v along 
the contracting CH. We shall show, analytically, that the 
divergence rate of both entities is exactly proportional 
to 1/r. Our method of calculation is non-perturbative, 
and is therefore valid also in the region of strong focusing; 
however, we shall use the perturbative results (applicable 
at the early section of the CH) to determine the two over- 
all coefficients characterizing the blue-shift divergence. 

The paper is organized as follows. In section || we de- 
scribe the physical model of the self-gravitating massless 
scalar field on a spherical charged black hole, and present 
the field equations. In Section III we carry out a leading- 
order perturbation analysis of r (and use pervious per- 
turbative results for <&) and calculate the w-derivatives 
of $ and r at the very early part of the CH (where the 
focusing effect is still negligible). Then, in Section |v| wo 
perform a fully nonlinear (and non-perturbative) calcu- 
lation of these ^-derivatives, which is valid everywhere 
along the contracting CH, up to the point of full focus- 
ing, where r = and the singularity becomes spacelike 
This nonlinear analysis leaves two coefficients unde- 
termined - one for each field - and we determine these 
two coefficients by matching the nonlinear results to the 
linear results applicable at the asymptotically-early part 
of the CH. Our results are in excellent agreement with 
the numerically-obtained results . 



II. PHYSICAL MODEL AND FIELD EQUATIONS 

We consider here the model of a spherically-symmetric 
charged black hole, nonlinearly perturbed by a self- 
gravitating, spherically symmetric, minimally-coupled, 
massless scalar field $ (the same model as that analyzed 
numerically in (|J^1). This model allows us to obtain non- 



1 



trivial radiative dynamics while retaining the simplicity 
of the spherical symmetry. 

We write the general spherically-symmetric line ele- 
ment in double-null coordinates, 



ds 2 = —f(u, v) du dv + r 2 (u, v) dO 2 



(1) 



where dfl 2 is the line element on the unit two-sphere. 
The energy-momentum tensor of a massless scalar field 
is 



T s = — I <T> <T> 

4?r " 



(2) 



The energy-momentum associated with the general 
spherically-symmetric free electromagnetic field is 
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where Q is the electric charge. 

For a spherically-symmetric scalar field, the Klein- 
Gordon equation reduces to 



(4) 



The Einstein field equations, = 8tt(T^ + Tffl, in- 
clude two evolution equations, 



p f,uf,v . j- f 
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(5) 



and two constraint equations: 

r,uu - (ln/) )U r> + r($„) 2 = 

r,w - (In /),„;> + r($ „) 2 = 0. 



(6) 

(7) 
(8) 



The form of the above line element and field equa- 
tions is invariant to a coordinate transformation of the 
form v — > — > u(it). In what follows u and u will 
denote generic, unspecified, double-null coordinates. Be- 
low we shall often use specific types of null coordinates 
for specific calculations, and in order to avoid confu- 
sion we shall assign a special symbol to each of these 
specific coordinates. Thus, we shall denote the stan- 
dard Eddington-Finkelstein null coordinates of Reissner- 
Nordstrom (RN) by u e and v e . We shall also use U and V 
to denote Kruskal-like coordinates, i.e. double-null coor- 
dinates which regularize the line element (0) at the inner 
horizon. In addition, in section IV we shall define two 



III. LINEAR REGIME 

Previous analytic and numerical studies have indicated 
that the geometry at (and near) the early part of the CH 
may well be described by the background metric func- 
tions of the static (or stationary) black-hole solution plus 
a small metric perturbation. This is found to be the situ- 
ation both in vacuum spinning black holes (analytically) 
H and in the present model of a spherical charged black 
hole (numerically) 0. Moreover, the perturbations be- 
come arbitrarily small as one approaches the asymptotic 
past of the CH. In the very early part of the CH, the 
perturbations are dominated by their linear part, and the 
singularity is well described by the linear metric pertur- 
bation. In the later part of the CH, however, nonlinear 
effects become exceedingly important, as demonstrated, 
e.g., by the contraction of the CH. 

Accordingly, we shall schematically divide the CH into 
two parts: 

1. The linear regime, i.e. the asymptotically-early 
part of the CH, where the metric perturbations 
(and the scalar field) are still very small, and a 
leading-order analysis is adequate; 

2. The nonlinear regime, i.e. the later part of the CH 
where the focusing (and possibly other nonlinear 
effects) become important. At the future end of 
the nonlinear regime the area of the CH shrinks to 
zero, and the singularity becomes spacelike. 

In this Section we shall consider the linear regime, and 
obtain expressions for the blue-shift factors, namely, the 
^-derivatives of <E> and r. The nonlinear regime will be 
the subject of Section Ev| . 

In the linear regime we may treat $ as a linear Klein- 
Gordon field over a fixed RN background. The evolution 
of such a field was analyzed in and more recently 
in [ ]13||l4| (using a different method). For a spherically- 
symmetric scalar field satisfying an inverse power-law 



$ = v e 



(EH) 



(9) 



at the event horizon (EH), the asymptotic behavior at 
the early part of the CH was found to be [ p^JT^ ] 



Bu K 



(CH), 



(10) 



where A and B are constants which only depend on the 
ratio Q/M, M being the black-hole mass. Since we are 
primarily interested here in the D-derivative of we only 
need the value of A, which was found in Refs. (l^Jl^l to 
be 



A =\ r -± r± + ^ 

2 r_ \r_ r + 



(11) 



other types of ingoing null coordinates, V r and V$. 



r± being the value of r at the outer and inner horizons 
of RN. 
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One finds that both at the EH and at the CH 

$ „ e oc v e ~ p , 
where p = n + 1, and 

$ CH /$ BH 



A. 



(12) 



(13) 



Here and below the arrow denotes the limit of large ad- 
vanced time (corresponding to v e — > oo). Note that the 
last relation is explicitly gauge invariant, so it holds for 
any type of ingoing null coordinate v, and not only for 
v = v e . 

Next, we consider the ^-derivative of r at the CH. In 
the pure RN geometry, r, e dies off exponentially (in v e ) 
at the CH. In the presence of the self-gravitating scalar 
field, however, r jVe decays as a power-law of v e (see be- 
low) . In the asymptotically-early portion of the CH (the 
"linear regime" ) which concerns us here, the effect of the 
scalar field is dominated by the second-order term (i.e., 
the term quadratic in derivatives of the scalar field), and 
higher-order corrections are negligible. We shall now cal- 
culate this leading-order term of r „ e . 

Viewing Eq. (^) as a linear first-order differential equa- 
tion for r „, we formally integrate it and obtain 



r, v (v) = -/(«) 



{V ' ] [$>')] 2 dv'. 



(14) 



Here [and also in Eq. ( Jl 5| ) below] the integration is done 
along lines of constant u, and we omit the dependence 
on u for brevity. From this exact expression we now 
extract the term quadratic in (derivatives of) $. Since 
$>, v 2 appears explicitly in the integrand, we simply need 
to replace / and r in the right-hand side by the corre- 
sponding unperturbed metric functions of RN, which we 



denote / RN and r RN : 



r, v (v) = -f R M 



[*„(«')] dv ' 



(15) 



Note that this equation is invariant to the choice of gauge 
for the coordinate v. We shall now evaluate the integral 
at the right-hand side, using the null coordinate v e . In 
this gauge, / RN decays exponentially at the CH: 



oc e 



(16) 



where k_ is the surface gravity of the inner horizon. On 
the other hand, $ >t , e decays as an inverse power of v e , 
and r RN approaches a nonzero constant, r_, at the CH. 
Therefore, since the relative change of (and r RN ) is 
exponentially slower than that of / RN , to the leading or- 
der in l/v e we can take <J? lK outside the integral, and 
substitute r RN = r_ [as well as Eq. ( |l6| ) for / RN ]. Doing 
so, we obtain (to the leading order in l/v e ) Q 



Finally, using Eq. (|l3|), we find 



In particular, we have in the linear regime 

T fY 11 * 



(17) 



(18) 



(19) 



One clarification should be made here concerning the 
precise meaning of the parameters r_ and k_ , and the 
coordinate v e , in the perturbed spacetime. (Originally 
these entities are only defined in the pure RN geometry.) 
We know that outside the black hole, both the scalar field 
and the metric perturbations decay at late time, and the 
geometry approaches that of RN. In particular, the mass 
function approaches a limiting value M. We thus define 
r_ and k_ according to the value of these parameters 
in the asymptotic RN geometry, i.e. according to their 
standard definition in terms of M and Q (with M be- 
ing the above late-time limit of the mass function; Note 
that the charge Q is a fixed parameter in our model). 
In a similar way, we also define the coordinate v e with 
respect to this late-time asymptotic RN geometry. More 
specifically, we may define v e according to the affine pa- 
rameter A along a line of constant r > 2M (or along the 
EH), by taking v e (X) to be the same function as in the 
pure RN geometry (with a mass parameter M defined as 
above). Note that once the entities M, r_, and v e 
were defined in the linear regime, their extension to the 
nonlinear regime is trivial. 

One might be puzzled by the relevance of the asymp- 
totic external mass parameter to the internal dynamics 
near the perturbed CH (and particularly to the definition 
of the inner-horizon parameters r_ and n~), especially 
when the divergence of the mass function at the CH is 
recalled. The resolution of this puzzle relays on the ba- 
sic features of the geometry inside perturbed charged (or 
spinning) black holes: On the one hand, the geometry 
is drastically different from that of RN (or Kerr), as ex- 
pressed by the divergence of curvature at the CH. On 
the other hand, the geometry is very similar to RN (or 
Kerr) in terms of the metric functions: The metric per- 
turbations are arbitrarily small at the asymptotic past of 
the CH. [Roughly speaking, the divergence of curvature 



* In the transition from Eq. (Tlil) to Eq. (|15[) we got rid of all 



terms of order higher than quadratic in $. Thus, in principle 
Eq. ( |l5| ) should include both the zero-order and the second- 
order parts of r. v . The zero-order term is represented by the 
(implicit) integration constant in Eq. (^|). This zero-order 
term is exponentially small, however, and is thus negligible 
compared to the quadratic term in Eq. (j(n\). 
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indicates the divergence of derivatives of the metric func- 
tions (with respect to the regular background coordinate 
V) at the CH.] This smallness of metric perturbations is 
the necessary basis for the entire perturbative approach: 
As it turns out, the perturbation analysis (when properly 
formulated) respects the smallness of the metric pertur- 
bations, and not the divergence of curvature. That is, 
the typical ratio of two successive terms in the nonlinear 
perturbation expansion is comparable to the small metric 
perturbations, and not to the diverging curvature (this is 
fortunate, because otherwise the perturbative approach 
would render useless). This was demonstrated analyt- 
ically for spinning black holes [0], and numerically for 
charged ones 0. 

In the above analysis of the linear regime (based on 
the perturbative approach), r_ and k_ appear as param- 
eters of the background RN geometry, and their definition 
should therefore be based on the asymptotic mass func- 
tion M. On the other hand, the divergence of the mass 
function (whose definition also involves the derivatives of 
r) at the perturbed CH merely reflects the divergence of 
r v there, due to the perturbation (which undergoes infi- 
nite blue-shift). Obviously, this divergence has no rele- 
vance to the background parameters r_ and k_ . 

IV. NONLINEAR REGIME 

We turn now to analyze the divergence rates of r, and 
along the nonlinear, strong-focusing, portion of the 
CH. Here, it will be insufficient to calculate the leading- 
order perturbations, so we must carry out a full nonlinear 
calculation. 

We shall base our calculation on two assumptions: 

1. For an appropriate choice of coordinates u, v, the 
line-element ([!]) is valid up to the singular CH, and 
both functions / and r are finite and nonvanishing 
along the singular CH. We shall denote such regular 
coordinates by U, V , and refer to them as Kruskal- 
like coordinates. (Of course, the choice of U and V 
is not unique.) We shall also set V — at the CH. 

2. There exists at least a single outgoing null geodesic, 
u = Mo, which intersects the CH and which satisfies 
the following two requirements: 

(a) Along u = uq, r and $ are monotonic func- 
tions of v in a neighborhood of the CH, 

(b) Along u = uo, both r v and d^/dr (i.e. 
$ v /r v ) diverge at the CH. 

The validity of assumption^ is strongly supported by the 
perturbative approach, at least in the early part of the 
CH. Moreover, recent numerical simulations Q confirm 
its validity in the entire CH up to the point of full focusing 



(where the singularity becomes spacelike). [] Assumption 
H is justified, because at least in the asymptotically-early 
part of the CH, Eqs. ( p^| , ^9| ) ensure the required mono- 
tonic behavior, and also imply dQ/dr oc v e p — > oo. In ad- 
dition, in the linear regime the standard ingoing Kruskal- 
like coordinate, V = e~ K - Ve , regularizes the line clement 
at the CH, and satisfies r v oc v e ~ 2p e K - Ve — » oo. We 
can thus take uo to be in this asymptotically-early sec- 
tion (in fact, the numerical simulations confirm that 
the asymptotic relations of assumption 2 hold everywhere 
along the CH). 

To analyze the evolution of r tV , we shall use the evo- 
lution equation (^), viewing it as a first-order ordinary 
equation for r „. Our goal is to integrate this equation 
in the ingoing direction, along the CH. This integration 
would be trivial if the last term on the right-hand side 
(which couples this equation to the other evolution equa- 
tion) were absent. Fortunately, on approaching the CH, 
this last term becomes arbitrarily small compared to the 
preceding one. For example, in a Kruskal-like V, the first 
term in the right-hand side diverges (at least at u = uq), 
whereas the second one is finite. (Note that although 
each of these terms depends on the gauge, their ratio 
is gauge-invariant.) This suggests that, when integrat- 
ing this equation along the CH, the last term could be 
dropped. In order to analyze this equation in a more 
systematic and elegant way, we define a new ingoing null 
coordinate V r in the neighborhood of the CH, by 

V r (v)=r(u = u ,v), (20) 

and reexpress Eq. (|5|) in terms of V r . To transform / = 
—2g uv from v to V r , we first calculate g uv : 

9 VVr = 9uvl (dV r /dV) = g uv / (dr/dV) Ug . (21) 

Since g uv is finite (assumption 1), and (dr /dV) Ug di- 
verges (assumption 2b), we find that g uv vanishes ev- 
erywhere along the CH, and so is g uV . Defining z(u) = 
(r Vr ) , Eq. (|^) now reduces to the trivial equation 
z. u = — (r u /r) z. Its general solution is 

z = C/r, (22) 

where C is an integration constant. Note that this exact 
equality holds everywhere along the CH. Calibrating C 
at u = Uo, we find 

^=7(^ = 7 (CH), (23) 



' If assumption yj were valid only in a portion of the CH, 
then the analysis below would nevertheless be applicable to 
this portion (provided that the outgoing null ray considered 
in assumption ^ intersects this portion). 
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where r$ is the r- value of the CH at u = uq. The first of 
these two equalities has an explicit gauge-invariant form, 
so we can immediately transform it to a generic gauge 
and write it as 



(r,.), 



'A 

r 



(24) 



(later we shall use this result for v = v e ). 

The analysis of the evolution of $ ;t , proceeds in a sim- 
ilar way. This time we use the KG equation ([|), viewed 
as an ordinary differential equation for <i> and integrate 
it along the CH. By virtue of assumption the second 
term in the parentheses in Eq. (||) is negligible at the CH 
(at least at u — uq) compared to the preceding one. To 
make an optimal use of this fact, we transform Eq. (^) 
from v to the new ingoing null coordinate 



V$(v) = $(u = tto, v), 



(25) 



defined in a neighborhood of the CH. The last term in 
the transformed equation is proportional to r v . But 



r„ (dV r /dV*)=r„ {dr/d<S>) l 



(26) 



At the CH, r lv = C/r and (dr/d$) Uo — > (assumption 
2b), so the last term in the transformed equation (^) 

vanishes. Defining y(u) = f$ v J , Eq. (|J) becomes 

* * ' CH 

y >u = —(r <u/r) y, whose general solution is 



y = K/r, 



(27) 



where K is an integration constant. Calibrating K at 
u = Uq , we find 



''() 



4> 



"0 



(CH), 



(28) 



and again, the first equality may be immediately trans- 
formed to a generic gauge: 



IE 

r 



(29) 



We shall now match the non- linear results (24) and 



(|29j) to the leading-order results at the linear regime. To 
that end, we take our reference outgoing ray u = uq 
to be in the asymptotically-early section of the CH. We 
can then use the results of the previous section [e.g. Eqs. 
( |l3| ) and dl8|) 1 for (r, v ) and {$, v ) Uo , and also substitute 
r = r_. Combining Eq. (|24| ) (with v = v e ) and Eq. (|29| ) 
with Eqs. (|l8[) and (|l3|), respectively, we obtain 



4> 



-.4 



and 



■ A 



(30) 



(31) 



These exact relations hold everywhere along the CH. 

More explicitly, for initial data $ = v e ~ n at the EH, 
the asymptotic behavior at the CH is (to leading order 
in l/v e ) 

„ 2 



r - A 

-n — Av e 
r 



-(n+l) 



■ A 



-2(n+l) 



(32) 



These results take an especially simple form when ex- 
pressed in terms of $ = r $ and r 2 : 

2 

^. Vc ^ -nr_Av e - {n+1) , (r 2 ) = -2n 



(33) 

(Note that to the leading order in l/v e , which concerns 
us here, the contribution of r tVe to *$> tVe is negligible.) 
That is, to the leading order in l/v e , the ^-derivatives of 
^ and r 2 at the CH are independent of r (and u). The 
translation of the above results from v e to any other type 
of ingoing null coordinate (e.g. V) is straightforward. 

The above results are verified numerically in Ref. [Q. 
The terms at the two sides of Eqs. ([30]) and (|3|) are 
evaluated numerically along an outgoing null ray that 
intersects the strong-focusing portion of the CH. The nu- 
merical results are in excellent agreement with the above 
analytic prediction. 

It would be an interesting challenge to try generalizing 
these results to the CH singularity of a generic spinning 
vacuum black hole. 
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